arXiv: 1508.01188vl [quant-ph] 5 Aug 2015 


Deterministic Quantum Computation With One Photonic Qubit 

M. Hor-Meyll, D. S. Tasca, and S. P. Walborn 
Instituto de Fisica, Universidade Federal do Rio de Janeiro, 

Caixa Postal 68528, Rio de Janeiro, RJ 2194F972, Brazil 

M. M. Santos 

Instituto de Fisica, Universidade Federal de Uberlandia, 

Caixa Postal 593, 38400-902, Uberlandia, MG, Brazil and 
Instituto de Ciencia, Engenharia e Tecnologia, Universidade Federal dos Vales do Jequitinhonha e Mucuri, 
Rua do Cruzeiro 01, Jardim Sao Paulo, 39803-371, Teofilo Otoni, Minas Gerais, Brazil 

E. I. Duzzioni 

Instituto de Fisica, Universidade Federal de Uberlandia, 

Caixa Postal 593, 38400-902, Uberlandia, MG, Brazil and 
Departamento de Fsica, Universidade Federal de Santa Catarina, 

Caixa Postal 4^0, 88040-900, Florianopolis, SC Brazil 

P. H. Souto RibeircQ 

Instituto de Fisica, Universidade Federal do Rio de Janeiro, 

Caixa Postal 68528, Rio de Janeiro, RJ 21941-972, Brazil and 
Departamento de Fsica, Universidade Federal de Santa Catarina, 

Caixa Postal 4^0, 88040-900, Florianopolis, SC Brazil 

We show that deterministic quantum computing with one qubit (DQC1) can be experimentally 
implemented with a spatial light modulator, using the polarization and the transverse spatial degrees 
of freedom of light. The scheme allows the computation of the trace of a high dimension matrix, 
being limited by the resolution of the modulator panel, and the technical imperfections. In order 
to illustrate the method, we compute the normalized trace of unitary matrices, and implement the 
Deutsch-Jozsa algorithm. The largest matrix that can be manipulated with our set-up is 1080 x 1920, 
which is able to represent a system with approximately 21 qubits. 


I. INTRODUCTION 

The field of quantum information has been strongly 
motivated by the demonstration that some quantum al¬ 
gorithms have an improved performance in comparison 
to their classical versions. More recently, the determinis¬ 
tic quantum computation with one quantum bit (DQC1) 
model was introduced with the aim of exploring the com¬ 
putational speed-up in high temperature ensemble quan¬ 
tum computation [I]. Although this model of compu¬ 
tation is not universal, it enables quantum speed up to 
solve certain problems, such as the Shor factorization al¬ 
gorithm [2], the measurement of the average fidelity de¬ 
cay of a quantum map [3], the trace calculus of an ar¬ 
bitrary unitary evolution [2., and the approximation of 
the Jones Polinomial [5]. The importance of the DQC1 
model of computing is that it requires little or no entan¬ 
glement between the qubits of the system || to evaluate 
the trace of a unitary operator, an operation that is not 
efficiently implemented by a classical computer [4]. The 
experimental implementation of this model of computa¬ 
tion has already been made in optical systems [7 and 
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We present an experimental scheme for the implemen¬ 
tation of the DQC1 protocol in an optical scenario. In 
this paradigmatic model, the normalized trace of a uni¬ 
tary matrix is computed by using the computational 
power of only one qubit in a pure state, and a collec¬ 
tion of qubits in a completely mixed state. The infor¬ 
mation about the normalized trace of the unitary matrix 
is transferred to the qubit state through conditional op¬ 
erations. For this purpose we use a phase-only Spatial 
Light Modulator (SLM), which performs polarization- 
controlled position-dependent phase shifts. Proper po¬ 
larization measurements return the result of the compu¬ 
tation. A similar approach was recently used to demon¬ 
strate the use of a SLM and polarized light to estimate 
integrals HOI. In addition to several possible applica¬ 
tions of the SLM El and its use as a quantum channel 
acting on the polarization [12], here we are interested in 
the implementation of quantum algorithms. Some quan¬ 
tum algorithms have already been implemented or simu¬ 
lated using the SLM, such as the Deutsch algorithm[13], 
Deutsch-Jozsa algorithm [HUS], the Grover algorithm 
m, and the quantum walk [16] . Moreover, there has 
been an implementation of the Deutsch algorithm for two 
qubits encoded in an optical system using an operation 
controlled by polarization El. but there an interferom¬ 
eter was needed to control the spatial mode. 

In order to illustrate the implementation of our 
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FIG. 1: DQC1 circuit. The action of the Hadammard gate 
on the state |0) gives the state |+) = (|0) + |l))/\/2, while the 
controlled unitary transformation U acts on the state p. 

method, some examples of trace calculation are presented 
for a few types of matrices and also the realization of 
the Deutsch-Jozsa algorithm. The oracle function is pre¬ 
pared in the SLM as a matrix of adjustable dimension, 
up to its resolution limits. For instance, for our High- 
Definition (HD) panel and using full resolution of 1080 x 
1920 pixels, we could implement an oracle function with 
approximately 2 21 inputs, which would correspond to an 
input composed of 21 qubits. We show that the perfor¬ 
mance of the method in the present realization is limited 
by polarization dephasing effects (12], and other sources 
of noise like fluctuations in the phase modulation and 
photon number statistics. 


the controlled operation, xi and yj are the coordinates 
of pixel i, j. Notice that in this description the state of 
the wavefront is maximally mixed. We will discuss later 
how this state can be obtained experimentally. 

Let us define an operator describing the action of the 
SLM on states of the form given by Eq. 0 as m- 

S=\H){H\®U+\V){V\®1, (3) 

where 

U='52e- UKxi ’«A\x i ,y j )(x i ,y j \, (4) 

id 

and cf)(xi,yj ) is a real function. 

The SLM used in our experiment only modulates 
the horizontal polarization component of the input light 
beam. The function <j)(xi,yj) is programmed to apply a 
phase between 0 and 27r in each SLM pixel located at 
(xi,yj). 

We use Eqs. to obtain the photonic state after 

incidence on the SLM: 

Pf = Sp i S' = S[\+){+\®Pt]& (5) 

= (l/2)(\H)(H\U Pt U^ + \V)(V\ Pt + 

+ \H)(V\U Pt l + \V)(H\lp t U'), 

and using the decompositions in Eqs. © and @ , we 
calculate the partial trace over the spatial degrees of free¬ 
dom: 


II. DQC1 USING A PHOTONIC QUBIT 

We start the discussion by explaining how we can use 
the SLM to implement the DQC1 circuit as sketched in 
Fig. [T] We use a light beam with a given transverse 
wavefront profile as the system p = p t , and its polariza¬ 
tion state as the control qubit. The initial photonic state 
can be described by: 


Ppo i = Tr t \p f ] = l(\H)(H\ + \V)(V\+ (6) 

+ \H)(V\Cj2*- iHxi ’ yi) + 

hi 

+ iv)mC£ 


Pi ~ |+)(+| ® pt , (1) 

where the polarization state |+) = (|0) + 11 ))/a/2 = 
(\H) + \V))/y/2, with | H) (|U)) representing the linear 
horizontal (vertical) polarization state and the state of 
the transverse wavefront is described by p t . 

Both the polarization and the transverse wavefront are 
properties of the same light beam. It is convenient to 
describe the wavefront in terms of a discrete basis: 


In terms of matrix representation the state can be writ¬ 
ten in the basis {\H),\V)} as: 

If 1 

Pp° l - 2 y CYsi •e +i ^ ( * <> ^ ) 1 ) - \ l ) 

This result shows that the information about the spa¬ 
tial modulation is transferred to the coherences of the 
polarization state, in terms of the average of the modu¬ 
lation distribution. Therefore, the expectation value ( a x ) 
of the Pauli operator cr x for this state gives: 


Pt = C^2\x i ,y j )(x i ,y j \, (2) 

id 

where each state of the basis describes the wavefront in a 
surface given by the pixel area of a SLM, and C is a nor¬ 
malization constant. The SLM will be used to implement 


(Px) =C^2 cos[<f>(xi,yj)\, (8) 

id 

which is the sum over the real parts of the modula¬ 
tion phases. We also have access to the sum over the 
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FIG. 2: Representation of the possible modulations for d 
2. White corresponds to a modulation phase 0, and blc 
corresponds to phase i r. 


imaginary parts of the modulation phases through the 
measurement of a y : 

K> = c Ys (9) 

We can interpret the modulation phases e~ l ^ Xuy ^ for 
each pixel as the diagonal elements of a matrix. There¬ 
fore, our scheme provides a method for calculating the 
normalized trace of this matrix through the measurement 
of ( a x ) and (a y ). 


A. Implementation of the Deutsch-Jozsa algorithm 

To illustrate further utility of the scheme, let us con¬ 
sider the implementation of the Deutsch-Jozsa algorithm. 
In this paradigmatic quantum algorithm, one wishes to 
test if an oracle function is constant or balanced m- 
The oracle function is implemented on the SLM. In the 
simplest case, also known as the Deutsch algorithm [20] . 
the SLM surface is divided in only two equal parts, where 
the upper part is given by yj > 0 and the lower part is 
given by yj <0. This corresponds to dimension d = 2, as 
illustrated in Fig. [2] For this algorithm, we use only two 
possible modulation phases for each half of the SLM: 0 
or 7r. For (t>(xi,yj) = 0, (a x ) = +1, and for <j){xi,yj) = 7r, 
(a x ) = —1. Therefore, for a constant function we have 
(cr x ) = ±1 and for a balanced function we have ( a x ) = 0. 

The algorithm is easily generalized (Deutsch-Jozsa al¬ 
gorithm) to any desired number of divisions of the SLM 
surface, up to the limit of a single pixel per division. In 
this case, we can have 1080 x 1920 = 2073600 cells that 
could be programmed individually with 0 or 7r. This is 
equivalent to testing a function with an input of approx¬ 
imately 21 qubits. The expected result for the Deutsch- 
Jozsa algorithm is the same as for 2 divisions, meaning 
that every kind of modulation map that modulates half of 


FIG. 3: Experimental set-up. He-Cd is a Helium-Cadmium 
laser, BBO is a nonlinear crystal, DM is a dichroic mirror, M 
represents mirrors, LI, L2 and L3 are lenses, A/4 is a quarter 
wave plate, A/2 is a half wave plate, BS is a 50:50 beam 
splitter, PBS is a polarizing beam splitter, SLM is the spatial 
light modulator, DET1 and DET2 are single photon counting 
modules, and CC is a coincidence detection circuit. 


the surface with 0 and half with i r implements a balanced 
function and should return zero as result. 


III. EXPERIMENT 

The experimental set-up is sketched in Fig. [3] A 
helium-cadmium (He-Cd) laser oscillating at 325nm is 
used to pump a type I BBO nonlinear crystal. It produces 
pairs of photons via parametric down-conversion, and 
we adjust the phase matching angle to obtain collinear 
twin beams at the degenerate wavelength of 650nm. The 
down-converted signal and idler beams are separated 
with a 50:50 beam splitter, and detected using lOnm 
bandwidth interference filters placed in front of the de¬ 
tectors. The idler beam is sent directly to a single pho¬ 
ton counter labeled DET1, and the signal beam is sent 
to a spatial light modulator (SLM) and polarization op¬ 
tics before detection by a single photon counter labeled 
DET2. Coincidence detection is used to post-select time 
correlated events signaling the arrival of a twin photon 
pair. In this way, the idler photon at DET1 heralds 
the presence of the signal photon. The signal photon 
propagates through LI, which is a lens implementing 
an optical Fourier transform mapping the far-held dis¬ 
tribution in the crystal plane onto the SLM. Using this 
lens system we avoid much of the free propagation ef¬ 
fects of the beams, minimizing unwanted diffraction and 
transverse spatial cross-correlations. The idler beam also 
propagates through LI before splitting in the beam split¬ 
ter. Its far-held distribution is mapped on an interme¬ 
diate plane, just like the signal but without an SLM. 
In our experiment, we use a rehective full HD 1080 x 
1920 pixels, phase-only SLM made by Holoeye Photon¬ 
ics. We are able to access each pixel individually and 
program any modulation phase ranging from 0 to 27r. 
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A specific function <p(xi^yj) characterizing the modula¬ 
tion distribution in the plane of the SLM is programmed, 
and the signal field wavefront acquires this phase condi¬ 
tioned on the polarization. The input beam is prepared 
in a linearly polarized state along the diagonal direction: 
|+) = (| H) + \V))/y/2, so that the horizontal compo¬ 
nent is modulated and the vertical is not. After interac¬ 
tion with the SLM, we measure the resulting polarization 
state given by Eq. ([7|, in two different basis, obtaining 
information about the modulation of the spatial profile. 
Lens L2 forms the image of the SLM surface in the de¬ 
tection plane. Lens L3 in the idler beam forms the image 
of the intermediate plane were the far field were mapped 
previously by LI. In this way, the spatial propagations 
of signal and idler beams are equivalent, except for the 
presence of the SLM in the signal. The polarization anal¬ 
ysis is made with a quarter wave plate A/4, followed by a 
half wave plate A/2 and a polarizing beam splitter (PBS). 
We register both the single photon and two-photon coin¬ 
cidence counts. 

Due to the spatial correlations, depending on how the 
idler photon is detected, the heralded signal photon can 
be prepared to have different spatial properties [21]. For 
a small idler detection area, the heralded signal beam 
becomes spatially coherent, which means that Eq. © is 
not suitable for describing its transverse wavefront. For 
a large detection area of the idler photon, the heralded 
signal beam is spatially incoherent, and is well described 
by Eq. ©• The use of twin photons and coincidence 
detection thus allows the preparation of states with con¬ 
trollable purity. 

It is important to characterize the spatial intensity dis¬ 
tribution of the light beam interacting with the SLM 
panel, as it may have some influence on the computa¬ 
tion. Fig. [4] shows the intensity profile of the signal light 
beam used, already triggered by the idler. This profile 
will be taken into account in the calculation of the trace 
of a matrix using our scheme. 



X T 


FIG. 4: Transverse profile of the beam on the SLM. (a) Sin¬ 
gle photon counting distribution (counts/s), (b) Coincidence 
counting rate (coincidences/s). The transverse coordinates 
Xi and Yj refer to the coordinates of detection cells contain¬ 
ing an array of (80 X 80) pixels. Cij is the single/coincidence 
counting rate in the cell located at (A/, Yj ) 


ing (1 x 1), (5 x 5), and (10 x 10) pixels. We can 
see that for all resolutions the measured value of ( a x ) is 
very close to zero, showing that the algorithm works well 
even when using the full resolution of the SLM. For the 
purpose of deciding if the oracle function is constant or 
balanced, a relatively high degree of uncertainty is toler¬ 
ated, since one is required to discriminate between 0 and 
±1. However, it is also clear that the error is smaller for 
larger sizes of the modulation cell, indicating the pres¬ 
ence of unwanted noise effects, like residual diffraction 
for instance. 


IV. RESULTS 

A. Implementation of the Deutsch-Jozsa Algorithm 

We begin by presenting the results obtained in the im¬ 
plementation of the Deutsch-Jozsa algorithm, as it is a 
special (and simpler) case of the general calculation of 
the normalized trace of a matrix. In Fig. [5] we show the 
surface modulation of the SLM and the measured values 
of ( a x ). For the constant functions the ideal value should 
be (<j x ) = +1 for 0 = 0 and (a x ) = —1 for </> = 7r, and we 
observe a good agreement with the measurements. For 
the balanced functions, we apply SLM masks with half of 
the cells modulated with 0 = 0 and half with 0 = tt dis¬ 
tributed randomly on the SLM surface. The randomness 
reduces the need of a precise knowledge of the spatial 
distribution of the light beam on the surface of the SLM. 
The results are shown in Fig. [5] for square cells contain- 


B. Computation of the Normalized Trace of a 
Matrix 

As we can see from the results in Eq. (|6|, we can 
implement the DQC1 model for a general function pro¬ 
grammed in the SLM, where we can define squares or 
other geometries composed by an arbitrary number of 
pixels and modulate each one with some arbitrary phase 
ranging from 0 to 27r. The measurement of ( a x ) gives the 
real part of the sum over all phases, and ( a y ) gives the 
imaginary part. Therefore, the system is able to perform 
the calculation of the trace of a normalized matrix im¬ 
printed in the SLM. In this case, it is not only a matter of 
determining if the function is balanced or constant, but 
rather of obtaining a result that is as precise as possible. 

In order to illustrate the method, we programmed the 
SLM panel with phases varying linearly along the y di¬ 
rection, according to: 
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Constant, 4> = 0 
<g x > — 0.92 


Constant, § = n 
<a x > = - 0.94 


Balanced, 1 pixel 
<a x > = 0.19 


Balanced, 5 pixels 
<g x > = - 0.04 


Balanced, 10 pixels 
<d x > = - 0.01 


FIG. 5: Images displayed on the SLM panel for constant 
and balanced functions and the respective measured values 
of {cr x ). The balanced functions are implemented through (1 
x 1), (5 X 5), and (10 x 10) square pixels cells randomly 
distributed on the panel. 






Vj) = 00 + ( 10 ) 


where N y = 1080 is the number of pixels along the y di¬ 
rection. Thus, the matrix element corresponding to posi¬ 
tion {xi^yj) is . We have considered four cases 

of linear functions of the type given in Eq. [lOj namely, 
with (0o,0/) = (37t/4, 57t/4), (7t,27t), (7t/2, 37t/2), and 
(7t/2,7t). The measured values (<r x ) exp and (a y ) eX p for 
the real and imaginary part of the trace of the corre¬ 
sponding matrices, respectively, are shown in Fig. [6] and 
summarized in Table I. 

We compared the experimental values ( (J x ) exp and 
(&y)exp with a theoretical prediction that takes into ac¬ 
count the non-uniform intensity distribution of the light 
on the SLM panel combined with the unavoidable de¬ 
phasing effect that is inherent to our Holoeye modulator, 
as described in detail in Ref. [120 The theoretical pre¬ 
diction for the real and imaginary parts of the trace of 
the matrix is well described by 


N X Ny 

{&x)theo Cij cos[<f>(xi,yj)], 

i =1 3= 1 
N x N y 

{p~y)theo — (1 - 2 p) EE Cij sinl^Xi, yj)], (11) 

i =1 3 = 1 


where p is a parameter which gives the degree of dephas¬ 
ing and is the measured intensity of the light incident 
on the SLM pixel at position {pc^yj). The overall de¬ 
phasing effect is to decrease the coherences in Eq. 0 
by a factor of (1 — 2 p). According to Ref. [12], the es¬ 
timated value for p, considering the specific modulator 
we are using, is 0.08 db 0.02. We determine Cij from the 
measurement of the transverse profile of the beam, shown 
in Fig. Qb). The resolution of our measurement of the 
beam transverse intensity is set by the size of a square 
cell composed by 80 x 80 pixels. Assuming a flat inten¬ 
sity distribution within the cell at position (A/, Yj), the 
intensity Q/ on a pixel at position {xi , yj ) inside the cell 
can be approximated by 


— 1 Cij 
' ij ~ N 80 x 80 ’ 


( 12 ) 


where Cu is the coincidence counts for the corresponding 
cell and A is a normalization factor given by 


N x /80 N y /80 

N = E E ( 13 ) 

i=i j= i 

so that 

N X Ny 

EE c *i = 1 - ( 14 ) 

i=1 3 = 1 


For the sake of completeness, the purely theoretical 
values for the real and imaginary parts of the normalized 
trace of matrix M (Tr n (M)), supposing the transverse 
field distribution to be flat in the entire SLM panel and 
disregarding the dephasing effect □a, given by 


Tr n (M) = Tr(M )/(N x .N y ), (15) 

are also shown in Table I. 

The error bars for the experimental values, shown in 
Table I, were computed considering the error in the phase 
modulation to be equal to the smallest modulation step, 
which is 5(f) ~ 27r/256, and considering the error in the 
intensity of the light in the SLM given by the uncer¬ 
tainty of the photon number distribution, considered to 
be Poissonian, so that Scij ~ y / G/- In general, we have 
a rather good agreement between experimental and the¬ 
oretical values within the technical limitations. For in¬ 
stance, they are related to the intrinsic discretization of 
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FIG. 6: Experimental and theoretical results of the compu¬ 
tation for the traces of the matrices defined by Eq. (10). 
Columns filled with diagonal(/) lines (blue) represent experi¬ 
mental results of measurements of (cr x ) (left side) and (a y ) 
(right side), and columns filled with anti-diagonal(\) lines 
(red) represent theory. The modulation ranges are displayed 
on the top. 


(00,0/) 

/ 3tt 5tt \ 

V 4 ’ 4 / 

(tt, 2tt) 

/ 7T 37T \ 

V 2 ’ 2 ) 

(f- 7r ) 

(&x)exp 

-0.811 

-0.039 

-0.646 

-0.579 


±0.005 

±0.008 

±0.006 

±0.007 

(&x) theo 

-0.773 

0.002 

-0.593 

-0.548 

Re[TV n (M)] 

-0.903 

-0.004 

-0.644 

-0.638 

(&y)exp 

0.007 

-0.628 

-0.034 

0.521 


±0.008 

±0.007 

±0.008 

±0.007 

(& 1 /) theo 

0.008 

-0.587 

-0.009 

0.545 

Im[Tr„(M)] 

0.012 

-0.637 

-0.003 

0.639 


TABLE I: Experimental and theoretical results for the com¬ 
putation of the traces of the matrices defined by Eq. (10). 


the modulator (here, we refer to the 256 levels of phase 
modulation), diffraction, inhomogeneities in the optical 
beam profile, and imperfections of the optical devices 
such as waveplates and polarizing beam splitters. Nev¬ 
ertheless, the results serve for a successful proof of prin¬ 
ciple. 


V. CONCLUSION 

In conclusion, we present an experiment in which 
the DQC1 quantum computation model is implemented 
using a polarization-controlled spatial light modulator 
(SLM) acting on the wavefront of single-photon fields. 
We illustrate the utility of the system by implementing 
the Deutsch-Jozsa algorithm on a system whose size is 
the equivalent of about 19 qubits using a matrix of 960 
x 960. The resolution of the SLM allows for the repre¬ 
sentation of the equivalent of up to about 21 qubits. We 
also show that a more general calculation is possible and 
experimentally compute, through polarization measure¬ 
ments, the normalized trace of a matrix whose diagonal 
elements are represented by modulation phases. 

A future path to improve our results concerns the use of 
alternative optical devices that could also implement po¬ 
larization controlled operations on the transverse spatial 
degrees of freedom of photons with better performance. 
One possibility could be the improvement of the SLM 
technology or the combined use of multiple SLMs. In 
this case we might be able to handle much larger matri¬ 
ces. 
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